Dipolar broadening of nuclear spin resonance under dynamical pumping 
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We study the polarisation dependence of the homogeneously broadened nuclear spin resonance in 
a crystal. We employ a combinatorial method to restrict the nuclear states to a fixed polarisation 
and show that the centre of the resonance is shifted linearly with the nuclear polarisation by up to 
the zero polarisation line width. The width shrinks from its maximum value at zero polarisation 
to zero at full polarisation. This suggests to use the line shape as a direct measure of nuclear 
polarisation reached under dynamical pumping. In the limit of single quantum of excitation above 
the fully ferromagnetic state, we provide an explicit solution to the problem of nuclear spin dynamics 
which links a bound on the fastest decay rate to the observable width of the resonance line. 

PACS numbers: 76.20.+q, 76.60.-k, 75.10.Jm 
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I. INTRODUCTION 

Recent interest in studying nuclear spins was fuelled 
[l|-[j| by the dephasing effect in GaAs quantum dots 
which is induced by fluctuations of the nuclear (Over- 
hauser) field and blocks prospective application of elec- 
tron spins for spintronics [9j and quantum computing 
[To| . Among different approaches to this problem, there 
is a new tool that combines radio frequency excitation 
and dynamical nuclear pumping techniques to control the 
Overhauser field and simultaneously to detect its state 



is reduced, relative to an unpolarised bath, as 
(iS)=3I(I + l)F 2 w(p), 



(2) 



11]. In particular, this method allows to access directly 
some of the local properties of the nuclear bath. 

At infinite temperature, a connection between the line 
width of the nuclear magnetic resonance and the dipole 
interaction between nuclei in a crystal was established 
in 12J. Later, it was found that the line shape changes 
at low spin temperatures [l3j]. A detailed understanding 
of the effect of a non-equilibrium polarisation, as occurs 
in dynamical pumping experiments, may give a way to 
measure directly the degree of polarisation reached un- 
der dynamical pumping and provide further insights in 
the effect of dipole interaction which plays an import- 
ant role in the intrinsic dynamics of the nuclear bath in 
unstrained structures. 

In present paper we address this problem by analysing 
the long wavelength spectrum of polarised nuclei using a 
cumulant expansion of the full line shape. We employ a 
combinatorial method to restrict the nuclear states to a 
fixed polarisation p = S z / {IN), where S z is the number 
of excitation quanta pumped to the nuclear bath, and 
find that the centre of the resonance shifts linearly with 
the polarisation, 



(u) = fiB + 37 Fp 



(1) 



where B is the external magnetic field, fj, is the nuclear 
magneton, / is spin of a nucleus, N is the total number of 
nuclei, and F is a maximum field induced by fully polar- 
ised nuclear bath. The line width at a finite polarisation 



where F 2 is a fluctuating field induced by completely un- 
polarised nuclear bath and w (p) is a bell shaped function, 
w (0) = 1 and w (±1) = 0; w (p) = 1 - p 2 for 1=1/2 
and the line width shrinks slightly faster for I > 1/2. 
The third cumulant is zero when p = but becomes 
comparable to the line width at a moderate polarisation 
signalling the build-up of a line shape asymmetry. 

One way to suppress the dephasing effect of the Over- 
hauser field is to polarise the nuclei above 99% [f|. So 
far the polarisation reached in dynamical nuclear pump- 
ing experiments is estimated indirectly via the hyperfine 
shift of the electron spin and, as reported, is still below 
the required threshold 0. The result of this paper sug- 
gests to use the shape of the nuclear magnetic resonance 
as a direct indicator of the degree of polarisation achieved 
in quantum dots. 

The vanishing of the broadening at large polarisations 
is caused by predominance of the exchange part of the 
dipole interaction. In the limit of single quantum of excit- 
ation above the ferromagnetic state the eigenmodes are 
plane waves with a sinusoidal dispersion, which allows 
to solve the problem of nuclear spin dynamics explicitly. 
We show that the observed line width of the nuclear res- 
onance in a single crystal of GaAs [13] gives a bound on 
the shortest polarisation decay times for lateral quantum 
dots as T\ > 1.25 s which is only an order of magnitude 
shorter than the T\ = 10 — 100 s measured in these dots 
2, 15] at p <^ 1. Note that, in this bound, we neglect 
quadrupole and hyperfine interactions of nuclear spins. 
Thus it is valid for unstrained structures without con- 
duction or localised electrons. 

The paper is organised as follows. In Section II we ap- 
proximate the model of dipole interaction between nuc- 
lear spins for a large external magnetic field. In Sec- 
tion III we analyse the shape the absorption line to ob- 
tain its polarisation dependence. Section IV contains the 
analysis of single quantum of excitation above the fully 
ferromagnetic state where we establish a bound on the 



2 



shortest poralisation decay time. In the appendices we 
give some details of the calculations in Section III. 



II. DIPOLE INTERACTION 

The Hamiltonian that describes N nuclear spins loc- 
ated at the sites of a regular lattice, their dipole interac- 
tion with the strength g, and an external magnetic field 
B is given by 

3 i<3 13 

(3) 

where fj, is the nuclear magneton, Ij,lf = ijii/j are the 
spin-/ operators, index j labels the sites on the lattice, 
r,j is the distance between two sites, and = v, L j 

Ay- 
is the unit vector that connects sites i and j. We con- 
sider a ci-dimensional crystal with N = L d nuclear spins 
and assume periodic boundary conditions, lj+L = Ij, to 
exclude edge states from the analysis. 

In a strong magnetic field the Hilbert space is parti- 
tioned by a large Zeeman energy into a set of subspaces 
labelled by the total z-projection of all spins, S z is an ei- 
genvalue of Ij ■ I n each subspace the many spin states 
that have the same Zeeman energy are split by a part of 
the dipole interaction that conserves S z , 



of the eigenstates |/) and |/') 



Hq — /.iB ' Sj 

3 



I -If - - J -^- 



i<3 



The matrix elements of the 
remainder, Hi = y\ . 3g 



4 (aa-ifa) 7ij iri?]/(8rl 



(4) 

S z non-conserving 

2 



connect 



subspaces 



with different Zeeman energies and can be treated 
perturbatively when the magnetic field is large. Here 
(dy , Pij , "/ij ) are coordinates of the vector e,j in a 
reference frame with z-axis parallel to the external 
B-field. 



III. ABSORPTION LINE SHAPE 

The absorption of electromagnetic radiation in a strong 
magnetic field is dominated by a flip of a single nuclear 
spin at an energy of /i-B, neglecting the dipole interaction. 
The dipole interaction of a single spin with other nuclear 
spins on the lattice makes the width of the resonance 
finite. Following [l2j we study moments of this line shape 
given by the transition energies Ej — Ef> between pairs 



fcx _ E//< (E f ,-E f ) k \(f'\I+\f)f 

1 ) ~ 



E//>iOT + i/>r 



(5) 



where I + = I~j i s the operator that corresponds to 
a long wave-length photon and k is the order of the mo- 
ment. To account for a finite polarisation, the initial 
states / are restricted to a subspace of fixed S z . In our 
study of the dipolar broadening we will use only the trun- 
cated part Ef\f) = Hq I/) of the dipole Hamiltonian H . 
A perturbative treatment of H\ produces extra peaks at 
2fj,B and 3fiB with small amplitudes ~ g/B, which are 
neglected. 

The sums in Eq. (JSJ) can be transformed to traces by 
means of general quantum mechanical identities. In the 
denominator, the sum over the eigenstates |/') gives a 
unit matrix, Ylf I/') (/'I = 1> m the subspace of S z + 1 
as the operator I + changes S z by 1 only. Then, the 
remaining sum over a complete set of the eigenstates is 
J2f (f\I~I + \.f) = Tr where the trace is restricted 

to a subspace of fixed S z . Representing the eigenenergies 
as Ef = (f\Ho\f) the numerator can be transformed in 
a similar way. Explicit expressions for the first three 
moments are 



(m) 



-Tr([g , /-]/+) 

Tr (/-/+) 
-Tr([ff ,J-][ff ,J+]) 
Tr (/"/+) 



(6) 



i3x _ Tr([H (h [H {h I-]}[H ,I+}) 



(m 3 ) = 



Tr (/-/+) 



where [A, B] = AB — BA is the operator commut- 
ator. Note, for calculating the absorption spectrum, that 
I + has to appear on the right of I~ in the traces. 



A. 1=1/2 case 

As a trace is invariant with respect to a change of basis 
we will calculate Eq. ([6]) in the basis of non-interacting 
spins that are quantised individually along the external 
magnetic field. First, let us consider the nuclei with spin 
1 = 1/2. In the basis of non interacting states, the de- 
nominator is a single sum of traces of individual spin 
operators Tr (jr/+) = £.Tr (ifi^). The full set 
of the basis states in a subspace of fixed S z can be di- 
vided in two groups: one has the i th spin in the state 
If = 1/2 and the other in the state If = -1/2. For 
each state form the first group the expectation value of 
If~Ii~ is and for the second group it is 1. Then, the 
trace of every single spin operator is independent of i, 



• C 



Tr(V* + ) 

factors C%-_\ and Cj^ give the numbers of states in 



N-l 



1 ■ Cff_i, where the binomial 
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each group. Here n is the number of 1/2-spins point- 
ing up with respect to the fully polarised state, i. e. 
S z — —N/2 + n. The sum over all N single spin traces 
gives the denominator as 



Tr (ri + ) 



AH 



n\ (N - 



1)!' 



(7) 



The numerator of the first moment k — 1 requires only 
one first order commutator, [Hq,I~] = — /J-B^^Li Ij ~ 
J2i<j (Sjfj — l) I'li- The trace of a two spin oper- 
ator, that is involved in Eq. ([6]), is evaluated analogously 
to the denominator. We calculate averages over all pos- 
sible two spin states accounting for the remaining N — 2 
spin states by binomial factors, 



Trinirz 



(1 - Sij) (N - 2)! (JV-2n-l) 



n! (AT- 



1)! 



(8) 



Substituting the commutator [77o,7~] into Eq. © and 
using the last expression we obtain the line shift, (v) = 
(m), that corresponds to the first moment 



3N-2n- 

(v) = uB - F 

w P 2 N - 1 



1 



(9) 



where F = gJ2j {^jfj ~ l) l r % is independent of the po- 
sition on the lattice i due to the periodic boundary con- 
ditions. Expanding the above expression in a 1/JV-series 
and identifying the polarisation as p = — 1 + 2n/N we 
obtain Eq. (p} for the case 7 = 1/2. 

The numerator of the second moment k — 2 requires 
an extra first order commutator [77 ,7 + ] = 7?^. 7 + + 
J2i<j rF^ (3zfj — l) Tjlf. The trace of a three spin op- 



erator, which is needed for Oi 
way, 



2 ), is obtained in a similar 



Tr(7/7r/^) = 



{I -5^) (1 - <5«/) (JV - 3)1 



4n! (N — n — 1)! 

[(2n -N+ l)(2n -N + 2)- 5 w 2n{2n - 2N + 3)] . 

(10) 

Substituting the product of the two commutators 
[H, I~\ [H, 7+] into Eq. ©, and using the above expres- 
sion, we obtain the second cumulant (v 2 ) = (m 2 ) — (to) 2 
that corresponds to the line width, 



<- 2 > 



F-, 



3n{2N-2n-3) 
2{N- l)(N -2) : 



(11) 



where F 2 = g 2 £\ (3^ - l) jr% r In the limit N > 1 
this expression leads to Eq. @ for the case 7 = 1/2 and 
w (p) — 1 — p 2 . 



The third moment contains an admixture of the first 
and the second moments when the polarisation is finite, 
p / 0. Hence, we evaluate the third irreducible moment 

using the same pro- 
In the limit N » 1, 



as 



(^ 3 ) = (m 3 ) - 3 {v 2 ) (v) - 
cedure as for the moments above, 
see Appendix A, we obtain 



9 



<- 3 > - § W 



1 3c~F 2 F-5F 3 )p(l-p 2 



(12) 



= <? 3 E ? 



H and Ft 



3c — 

This 



where 7*3 

9 3 E« (H ~ l) (375 ~ 1) (3rg - 1) / (r&r&rg 
result indicates that the line shape is asymmetric when 
the nuclear spins are polarised. The asymmetry be- 
comes pronounced when the skewness (^ 3 ) / (V 2 ) 3 ^ 2 = 
y/3/F$ (F 3c - F 2 F - 5F 3 )p/^l-p 2 is of the order of 
one. For example, in a cubic crystal with external mag- 
netic field B || [001] this point is reached at a mod- 
erate polarisation of p = 0.32. 



B. 



I > 1/2 case 



For nuclei with 7 > 1/2, the extension of the combin- 
atorial calculation presents an extra difficulty due to the 
larger number of single spin states. We start from the 
denominator of Eq. © which involves a trace over only 
one single spin operator, Tr (jr JT ) . Summing over 27+1 
single spin states, instead of just If — ±1/2, modifies Eq. 
to 



21 



Tr (7-7+) = (k + 1) (27 - k) 



(13) 



k=0 



IS 



where lA » = E^ 1 ' (-1)* ^ C^^f^ 
the number of linearly independent eigenstates of the op- 
erator E • 7| with S z = —NI + n, constructed out of A" 
spins 7, see Appendix B. This provides a generalisation 
of the binomial factor for 7 > 1/2; it can be check that 



A 



N 



The traces in Eqs. ([H [T0| of the two and the tree single 
spin operators which are needed for the first two moments 
are modified in the same way, using jA% instead of the 
binomial factor, 



2/ 



k 2 

fci,fc 2 =0 

x (-7 + fc 1 )(fc 2 + l)(27-fc 2 ) 



(14) 



4 



27 

Tr (1*1-1+11) = (1 - Sa) (1 - S u ,) £ iK- k \- k2 (^1 + 1) (2/ - ^1) {-I + fc 2 ) 2 

fei,fe 2 =0 
27 

+ (1 - 8 jV ) (1 - (1 - 8 iV ) ^ (fci + 1) (21 - fci) (-/ + fc 2 ) (-/ + h) (1 - 5 W ) (15) 

fcl ,fc 2 ,fc3=0 



We substitute Eqs. (|13I14I15|) into Eq. © and evaluate 
the cumulants numerically. The first cumulant coincides 
with Eq. ([1} in the limit of a large crystal N 3> 1, while 
the second cumulant gives the Eq. @ with w (p) evalu- 
ated in Appendix C. 

In a crystal with two species of nuclear spins the di- 
polar broadening of the like spins has an extra contribu- 
tion from the unlike spins. We assume that the corres- 
ponding magnetons, [i a and are significantly different, 
I Ma - W>| B > g a> , g b , g ab where g a ,,9b, 9ab are the dipolar 
coupling constants between the like and the unlike nuclei. 
Then the off-resonant interaction between the unlike nuc- 
lei makes the exchange term of Eq. (TJ| a perturbation, 
which we will neglect. For the unlike spins we repeat 
the same calculation as for the like spins keeping only 
the Ising term of the dipolar interaction between the two 
sublatices. We find that the first two cumulants are sums 
of two contributions, 

(v a ) = » a B + (l + 2S ac )I c F ac p Cl (16) 

c—a,b 

^ = Y / 1 -±^I c (I c + l)F^w(p c ), (17) 



where F a 



9aa Jji 



l) is a sum over sublat- 
— 1) is over sublattice 



tice a and F ab &V.^(.. 

b; both sums are evaluated in the same reference frame. 
The polarisations of each species, p a and p b , are kept as 
independent parameters. 

The line shift Eq. ((T|) can increase the line width meas- 
ured via the hyperfine shift of the electron Zeeman line 
In this technique all nuclei are polarised first, then 
they are depolarised at one frequency after another and 
shifts of the electron line are measured. But in the depol- 
arisation process the line position changes that can add 
to the measured line width. This addition is anisotropic 
similarly to the zero polarisation line width in Eq. @. 
For Zinc Blende crystal of a III-V semiconductor with 
I„. = If, = 3/2, using Eqs. ([TBI fl7|) the maximum shift 

=i " (fa)p=o) /y/Jffi is when B || [001] and is 
0.8 when B | [111] for one sublattice [l6[ assuming that 
the other, off-resonance, sublattice stays polarised. 



(< 



IV. SINGLE QUANTUM OF EXCITATION 
ABOVE THE FERROMAGNETIC STATE 

Now we consider the situation where nuclei are close to 
fully polarised. The line width Eq. © decreases at these 
large polarisations because the exchange term in Eq. Q 
becomes dominant. The eigenmodes of the exchange 
term, the magnons, are plane waves, so only a few, with 
long wavelengths, overlap significantly with the radio fre- 
quency photon in Eq. (|5|) . In the limit of single quantum 
of excitation above the fully ferromagnetic state, corres- 
ponding to the subspace of S z = —IN+1, the eigenmodes 
of full Hq are plane waves. As an example, for a cu- 
bic crystal they are |q) = J2j ex P (Q ' r ij) liO /V2IN 
where the momentum is quantised by the periodic bound- 
ary conditions q = 27rQ/ (L + l) d and < Qi < L + 1 
are integers. The dispersion is obtained directly using 
the Schroedinger equation, H Q \q) = (E^ + E\) |q), as 

= -^Ej:,. 3 >0 COS (q- r ij)( 3 7 l 2 J -l)/4- Here ^1 

is a momentum independent part. The transition mat- 
rix element in Eq. §5§ for the initial state [/) from the 
subspace S z — —IN is the largest for the smallest mo- 
mentum q = 2n/ (L + l) d and is suppressed for all other 
momenta giving the zero line width in the limit N 3> 1. 

In the subspace of S z = —IN + 1 the simple form of 
the eigenmodes permits explicit calculation of the spin 
dynamics. The evolution of a locally created polarisation 
|0> 

spreading profile at later times 



/"l/2 |4) /y/2I, due to the action of H , gives a 



(t\i;\t) = -i- 



E 



iq.(r— |)-iE q i 



N 



(18) 



For a large magnetic field, corrections to this average due 
to the off-resonant part of the dipole interaction Hi will 
involve the two and the three excitation subspaces. In 
the specific case of d = 1 and only nearest neighbour 
interactions, the Bethe-Ansatz approach [19( allows to 
assess these corrections, which turn out to be small, ~ 

g/B. 

The width of the probability profile Eq. (fT8)) grows in 
space linearly with time r ~ tg/A 2 [3], where g is the di- 
polar coupling g and A is the lattice parameter. At smal- 
ler polarisations \p\ < 1 the Ising interaction leads to a 
magnon-magnon scattering which slows the free magnon 



5 



propagation. This sets a bound on the shortest nuclear 
polarisation decay time for a quantum dot of a radius tq 
as Ti = r A 2 /g. The dipolar coupling constant can be 
extracted from the observed width of the nuclear mag- 
netic resonance. In a single crystal of GaAs the full width 
at half maximum is 5 kHz [14|. Using Eq. (TTT| for zero 
polarisation and assuming gn^a^ = g As As = ffGaAs = 9 
we extract g = 12 Hz • nm 3 [16(. For a lateral quantum 
dot of the radius ro — 50 nm the bound is thus T\ > 1.25 
s. 
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Appendix A: Third moment for 1—1/2 
Here we evaluate the third cumulant as 



CONCLUSIONS 



(Al) 



In conclusion, we have analysed the effect of a non- 
equilibirum polarisation on the nuclear magnetic reson- 
ance and have found a direct relation between the line 
shape and polarisation. This result suggests to use the 
line shape to measure directly the degree of polarisation 
achieved in dynamical nuclear pumping experiments. 



The first and the second cumulants are evaluated in Eqs. 
(9,11)- The third moment in Eq. (6) involves a product 
of two commutators 

[H ,I+] = *£'/ + £ £3 (34 -1)////, (A2) 

j i<j v 



J 



[H ,[H ,r]]= £ -f-3(3 7 5-l)^-(3 7 ^-l)Q 



2^ 



l 1 ^ 



r,rr 
33 * 



(A3) 



We evaluate the trace of the operator parts of each term in the product separately. Keep in mind that the terms 
with i = j or j = j' do not enter into the sums thus we do not consider them. For the following terms we add up 
contributions from all possible configurations of three different spins 



-Tr I + I~I~,IZI Z „ = 



(A4) 



The above term is zero as we do not account for j = i and j = j' terms. 



-Tr (jflj Ij'lplf'?) — ^ij'^33"J ipN-3 ~ + (^TV-3 ~ @N-3 + ^N-z)) ' 



(A5) 



For the following terms we add up contributions from all possible configurations of four different spins 



Tr (ip*lrl+l^ 



n-2 
JV-3 



/in 
°iV-3 



c 



■n— 3 
AT-4 ' 



3C 



'71— 1 

N-i 



til" 
2 3 

1 <; J s~in— 3 1 /~in q/-m— 1 1 q/^ro— 2 /~m 

-\-0j'i" \ — + u JV-4 — ,JO JV-4 ' ,30 JV-4 °JV-3 



■ ( fin— 3 /in 
+<!>.,-,:» <%<{ - C 



3C 



A r_ 4 — JO JV _ 4 j- 



n-2 
N-3 



/~~tn—3 
L 'JV-4 



u iV-3 



AT-2 



/~in— 1 
u iV-2 



c 



Oy^m— 1 

0< ^Ar_4 



+ C"I 3 + 2C 



2C 



3C^_ 2 4 } 
t- C^_ 2 3 } 



1 

Af-3 



3C n-2 



N-A 



N 



N-3 



^N-4 



2 \ 
N-41 

2C 



1 

W-3 



(A6) 



q/"m— 2 \ 



6 



2Tr 



(9 



rz r- r+ rz 
2 J j V^' 



Sj'iSii" {CJJr_4 — C^r_ 



2 2 



-zo Ar _ 3 -t- oo A ,_ 4 - 



sin— 1 /^m r )C in 



AT-2 



+4 



f (~m—2 



3 

JV-3 " 
n-2 
3 



{c n N 
'{c n N -J 

>8u»{C 



N-3 

Gn 
AT-3 

-c 



AT-3 



os~tn—2 \ 

/-m—Z , 
U JV-4 "+ 

/^in— 3 
~ °JV-4 



JV-4 



3C 



JV-4 



C 



N-4 



3C% 



N- 



in 

o^-m— 1 
0o AT-4 



n— 3 
AT-4 ' 



c 



n - 



c 



n 



c 



N- 



3 ~ ^W-4 ^W-4 ; 

i o/~m — 1 q/~^ — 2 \ 
4 T ' 3L 'Ar-4 — • 3 °JV-4/ 



3C^ 



'AT-3 

nrtn—2 \ 
•J^JV-4/ 

<ys~in— 2 \ 

3C 



"\°iV-4 — °JV-4 ' 0L/ A'-4 ' 



6Lj N-4i 



/~in— 3 
°Af-4 



c 



Af-4 ' 



o/~in— 1 
0l - / Ar_4 



3C^_ 2 4 } 



(A7) 



Then, we substitute the combinatorial factors explicitly and divide them by the denominator Eq. (7). Keeping only 
the leading terms for 1/N -C 1 we get four contributions to the third moment, 



2 Tr (/-/+) 



n (-N + 2n) 2n{N - n) 

■ + OH' 



N 3 



N 3 



(A8) 



Tr 



(9 



3 « V' 



Tr (J-/H 



''7 

8 



4n(2n 2 - 3nJV + TV 2 ) 



(-iV + 2n} : 
N 4 



iV 4 

2n(-NN- 



4n(2n 2 - 3nN + N 2 ) 



2n) 2 



N 4 



8ji" Sij 



N 

8n 2 (n — N) 



(A9) 



N 4 



Tr 



/ TZ TZ J- T+ TZ 

Tr (/-/+) 



J 33 

4 



4n (2n 2 - 3niV + TV 2 ) An (2n 2 - 3nJV + N 2 

Sj'iSw jrg Sj'i 



N 4 



(A10) 



An (2n 2 - 3nN + N 2 ) An (2n 2 - 3nN + N 2 ) 

Su" ^4 S yi " — A 

An (2n 2 — 3iiN + N 2 ) (-N + 2nf 

+0jH0jH" h Sj'i"6ii" 

(-N + 2nf (~N + 2n) 

H — o-j>i»dii»o. 



N 4 



3n 



iV 4 



N 4 



Substituting these expressions in to the sums Eqs. 
(|A2IA3j) and identifying p = -1 + 2n/N we obtain Eq. 
(12). 



Appendix B: Derivation of iA n for arbitrary I 

Here we derive the expression for the combinatorial 
coefficient for arbitrary / which uses only one sum over 
the binomial factors. We start from the generating func- 
tion for a single spin 

f(x) = l + x + ---+x 21 = i—^ . (Bl) 



Then the combinatorial factor is given by the coefficient 
in front of the n th power of x in 

JV oo 

f N {x) = ^^(-l^^+^^C^ 1 ^^ 

oo 

= J2 (-l) fcl C£<" 1+fe2 z (2/+1)fel+fc2 (B2) 

fei,fe 2 =0 

where (l + a:)^ - Ek=o C k^ k and (1 - X y N = 
"l2T=o C^~ 1+k x k were used. Here the coefficient is 

n/(2/+l) 

AN _ ( 1 \ k r<N r ,N-l+n-(2I+l)k 

I A n - 2^ ° fe °n-(2/+l)fe ■ 

fe=0 
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Appendix C: Numerical evaluation of the first two 
moments for 7 > 1/2 

For a nuclear spin I > 1/2 the cumulants of the line 
shape are expressed in terms of sums over the binomial 



coefficients. We analyse these sums numerically. Sub- 
stituting Eqs. (13,14) into Eq. (6), we obtain the first 
moment, 



(v) -hb _ Eli M = iK-L-k. (- 1 + fc i) (** + 1) (21 k 2 ) 



iN-2 



F 



(CI) 



This expression coincides with Eq. (1) when N 1, see Figured] 

The polarisation dependence, w (p) = (v 2 ) / (31 (I + 1) F2), of the second cumulant is obtained by substituting Eqs. 
(12,13) into Eq. (6), 



w(p) 



^V 2/ tA 

° Z^fci,fe 2 ,fe 3 =0 IIi n-k 1 -k 2 -k 3 



tN-3 



(fcl + 1) (21 - fci) (-/ + fc 2 ) (-/ + fc 3 ) 



ElLoiK-k(k + l)(2I-k) 



3 E^, fe2 , fe3 =o i<-t- fe - fe3 + 1) ( 2/ - fc (- 1 + M (-/ + h) 



lJV-3 



+ 



EfeioMn-t(fc + l)(2/-fc) 

I 



(C2) 



■ Eq. (1), S=3/2 
Eq. (CI), 5=3/2 

• Eq. (1), S=7/2 

- - Eq. (CI), 5=7/2 




Figure 1: The first moment of the line shape. The full and 
dashed lines are linear function from Eq. (1) for I — 3/2 and 
I = 5/2. The squares and ellipses are Eq. ()C1I) numerically 
evaluated using for I = 3/2 and I = 5/2; N = 500. 



— 1-f 

• Eq. (C2), S=l/2 

-.- Eq. (C2), S=3/2 

-♦- Eq. (C2), S=9/2 




Figure 2: The second moment of the line shape. The full line 
is w (p) — l—p 2 . The ellipses, squares, and diamonds are Eq. 
(fC2| numerically evaluated for 1 = 1/2, I = 3/2 and / = 9/2; 
N = 500. 



In the limit 1/N <C 1 this expression coincides with 
w(p) = 1 - p 2 for I = 1/2. When I > 1/2 the nar- 
rowing is slightly faster than l—p 2 , see comparison in 
Figure. [5J 
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